Orthorhombic fulleride (CH 3 NH 2 )K 3 C 60 close to Mott-Hubbard instability: Ab initio 

study 



(N 

o 

(N 

C 
3 



-4— > 

CO 

a 

C 

o 
o 



> 
in 
m 
m 
(N 

O 



Anton Potocnik/'B Nicola Manini, 2 ' 3 ' 4 '0 Matej Komelj, 1 Erio Tosatti, 3 ' 4 ' 5 ^ and Denis Arcon 1 ' 6 

1 Condensed Matter Physics Department, J ozef Stefan Institute, Jarnova cesta 39, SI- 1000 Ljubljana, Slovenia. 
2 Dipartimento di Fisica, Universitd degli Studi di Milano, via Celoria 16, 20133 Milano, Italy. 
3 Scuola Internazionale Superiore di Studi Avanzati, Via Bonomea 265, 34136 Trieste, Italy. 
4 CNR-IOM Democritos, Via Bonomea 265, 34136 Trieste, Italy. 
5 International Centre for Theoretical Physics (ICTP), Strada Costiera 11, 34151 Trieste, Italy. 
6 Faculty of mathematics and physics, University of Ljubljana, Jadranska 19, SI-1000 Ljubljana, Slovenia. 

(Dated: March 8, 2013) 

We study the electronic structure and magnetic interactions in methylamine-intercalated or- 
thorhombic alkali-doped fullerene (CH 3 NH 2 )K 3 C 60 within the density functional theory. As in the 
simpler ammonia intercalated compound (NH 3 )K 3 C 60 , the orthorhombic crystal-field anisotropy A 
lifts the tiu triple degeneracy at the Y point and drives the system deep into the Mott-insulating 
phase. However, the computed A and conduction electron bandwidth W cannot alone account for 
the abnormally low experimental Neel temperature, Tn = 11 K of the methylamine compound, com- 
pared to the much higher value Tn = 40 K of the ammonia one. Significant interactions between 
CH 3 NH 2 and Cgg are responsible for the stabilization of particular pseudo-Jahn- Teller fullerene- 
cage distortions and the ensuing low-spin 5=1/2 state. These interactions also seem to affect 
the magnetic properties, as interfullerene exchange interactions depend on the relative orientation 
of pseudo-Jahn- Teller distortions of neighboring Cgg molecules. For the ferro-orientational order of 
CH 3 NH 9 -K + groups we find an apparent reduced dimensionality in magnetic exchange interactions, 
which may explain the suppressed Neel temperature. The disorder in exchange interactions caused 
by orientational disorder of CH 3 NH 2 -K + groups could further contribute to this suppression. 

PACS numbers: 71.20.Tx, 74.70.Wz, 71.30.+h 



I. INTRODUCTION 



The alkali-doped cubic fullerides (A 3 - X A' X C 60 , A, A' = 
Na, K, Rb, Cs) are prominent members of a rapidly grow- 
ing family of 7r-electron organic superconductors^ with 
superconducting transition temperature, T c , as high as 
38 In the early days of fullerene research the ob- 

served monotonic increase of T c as a function of unit cell 
volume was attributed to the increased density of states 
at the Fermi level, within a standard Bardeen-Cooper- 
Schrieffer (BCS) theory^ Due to the three-fold degener- 
acy of the ti u -derived bands at the T point, preserved in 
cubic structures, and to the well-established vibron cou- 
pling, the Jahn- Teller (JT) effect is believed to partici- 
pate actively in the superconducting pairing mechanism^ 
However, since the conduction electron's bandwidth (W) 
and the on-site Coulomb repulsion energy (U) are com- 
parable it was also argued that fullerides must be very 
close to a Mott-Hubbard transition^ Phase transitions 
from a superconducting to a Mott insulating state were 
indeed demonstrated long ago upon intercalation-induced 
lattice expansion in compounds such as (NH 3 )K 3 C 6C p 
where, similarly to Li 3 (NH 3 ) 6 C 60 ^ expansion increases 
the repulsion/bandwidth ratio U/W raising the impor- 
tance of correlations. These notions recently received 
strong support with the discovery of Cs 3 C 60 , a cubic ful- 
leride with the largest unit cell and low-spin (S = 1/2) 
antiferromagnetic Mott-insulating ground state at ambi- 
ent pressure conditions^ - — and where high-tcmpcrature 
superconductivity is revived under hydrostatic pressure 



with no structural symmetry change. Strikingly, under 
pressure here T c displays a dome-like unit-cell volume 
dependence typical of unconventional superconductors^ 
frankly inexplicable within BCS, Migdal-Eliashberg and 
related weakly-interacting models. Conversely, precisely 
such a nonmonotonic behavior of T c had been predicted 
by dynamical mean field theory (DMFT) for the un- 
conventional superconducting phase bordering a Mott- 
insulating state at large unit cell volumes in a simple 
three-band Hubbard model incorporating, besides a JT 
electron-vibron coupling, strong on-site electron correla- 
tions, caused both by the Coulomb repulsion U and by 
Hund's rule exchange, Jh^ That success not only high- 
lights the unconventional interplay of strong correlations 
and phonons in these systems: it also brings expanded 
alkali-doped fullcrenes at the forefront of the strongly- 
correlated electron systems and superconductors, where 
good understanding of one model system may help shed 
more light on the entire highly controversial field. 

In search for higher transition temperatures, ear- 
lier approaches toward expanded unit cells had been 
by means of co- intercalation of inert ammonia (NH 3 ), 
e.g. NH 3 K 3 . x Rb x C 60 (x = 0, 1, 2, 3)^ and in 
Li 3 (NH 3 ) 6 C 60 ^ or of methylamine (MA) molecules, in 
(CH 3 NH 2 )K 3 C 60 (MAK 3 C 60 )3 These molecules oper- 
ate as spacers between C 60 molecules and are believed 
to have negligible direct influence on the electronically 
active t\ n molecular orbitals of C 60 3_ anions. However, 
their presence breaks the original cubic symmetry result- 
ing in orthorhombic crystal structures and effectively lift- 
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ing the three-fold degeneracy of the ii u -derived bands at 
the r point. This splitting reduces the critical value for 
U/W c ,— immediately pushing non-cubic fullerides over 
the metal-insulator transition (MIT) boundaryiSi We 
stress that these anisotropic compounds would still be 
metals, presumably also BCS superconductors, if strong 
on-site electron correlations were not present^ They are 
thus important model systems for investigations of cor- 
relation effects in metals where the splitting of the origi- 
nally degenerate narrow bands can be tuned through the 
anisotropy field. 

At low temperatures NH 3 K 3 x Rb x C 60 compounds or- 
der to an antifcrromagnctic insulating state with Neel 
temperature, Tn, ranging from 40 K to the maximum of 
80 K for x = and x = 2, respectively.-^- MAK 3 C 60 , 
on the other hand, orders to an antiferromagnetic state 
at remarkably lower Tn = 11 Assuming that the 

electronic properties depend mainly on the unit-cell vol- 
ume one would anticipate Tn to be ~80 K, i.e. sim- 
ilar to that of NH 3 KRb 2 C 60 . Such profoundly differ- 
ent magnetic response to different co-intercalands is un- 
expected and implies that in addition to direct inter- 
fullerene electronic overlap also other degrees of freedom 
play a role in non-cubic fullerides. We first note that 
the c/a lattice parameter ratio is smaller in MAK 3 C 60 
than in NH 3 K 3 _ x Rb x C 60 , where it is closer to unit y 16 i 18 
Therefore, the difference between the two orthorhom- 
bic compounds may arise from the different orthorhom- 
bic anisotropy field, A. In other words, properties 
of MAK 3 C 60 and NH 3 K 3 _ x Rb x C 60 compounds are fine- 
tuned by the closeness to the Mott-Hubbard boundary 
on the (U, A) Manini-Santoro-Dal Corso-Tosatti (MSDT) 
phase diagram for fullerides ^ While it is well established 
that NH 3 K 3 C 60 lays very close to this boundary, the pre- 
cise position of MAK 3 C 60 is yet to be determined. 

In cubic fullerides the high crystal symmetry prevents 
ordering of the ti u molecular orbitals, i.e. it is believed 
that cubic A^Cqq are in an orbitally liquid stated at least 
above the antiferromagnetic Neel temperature. In the 
non-cubic fullerides, however, one expects that crystal 
field will select a particular Cg cage distortion and thus 
would also affect the inter-fulleridc hopping integrals. 
Since MA molecule is rather large compared to avail- 
able octahedral space in the fulleride structure a strong 
interaction - presumably even the formation of a weak 
hydrogen bond - between MA-K+ group and the near- 
est C 60 was proposed J£ MAK 3 C 60 thus offers a unique 
opportunity to investigate the role of symmetry breaking 
on the JT effect in strongly-correlated electron systems. 
Stimulated by these open issues we decided to carry out 
a systematic ab initio study of MAK 3 C 60 by the density 
functional theory (DFT) in the local density approxima- 
tion (LDA) , where even if strong correlations are treated 
only at the mean-field level, the detailed chemical bond- 
ing and crystal-field strengths can be assessed. We show 
that the crystal-field anisotropy of MAK 3 C 60 is larger 
compared to NH 3 K 3 C 60 , thus placing it deeper in the in- 
sulating region of the phase diagram. The whole electron- 



hopping structure is found to be significantly affected by 
the MA insertion. We find evidences for the presence of 
strong Cg pseudo-JT effect that may lead to a reduced 
dimensionality in magnetic exchange interactions and ex- 
plain suppressed Tn- 



II. METHODS 

We execute DFT calculations using the Quantum 
Espresso software package (pwscf program)^ Ultrasoft 
pscudo potentials appropriate for the Perdew-Zunger 
exchange-correlation (LDA) are used. The pseudo Bloch 
functions arc expanded over plane waves with an en- 
ergy cutoff of 50 Ry on a 5 3 Monkhorst-Pack A;-space 
mesh, such that the total energy is within 250 mcV of its 
converged value, while the conduction bands arc within 
0.1 meV of their converged value. The Fermi surface is 
smeared with a "temperature" parameter of 2 meV. The 
density of states (DOS) is evaluated with the tetrahe- 
dron method, sampling A;-space with a uniform 8 3 mesh. 
For the relaxation of the atomic positions, we adopt a 
damp (quick-min Verlet) procedure on a 2 3 fc-space mesh 
and plane-waves' energy cutoff of 30 Ry. The maximally- 
localized Wannier orbitals are obtained using the Wan- 
nier90 packaged The three Wannier orbitals are com- 
puted in the ti u energy window, using 5 3 Appoints on 
Monkhorst-Pack grid with Bloch phases as initial pro- 
jections. All three Wannier orbitals are positioned on 
the central C 60 molecule and have a similar spread of 
\[Q ~ 0.4 nm. The projected DOS (PDOS) on the 
three Wannier orbitals were evaluated using the Gaus- 
sian smearing method with a "temperature" parameter 
of 8 meV. The use of two different methods for evalua- 
tion of DOS and PDOS gives rise to a slight mismatch 
between the two. The mismatch does not affect the con- 
clusions of this work. 



III. RESULTS AND DISCUSSION 

A. Phase diagram 

The MAK 3 C 60 compound grows in the face-centered 
orthorhombic (space group Fmmm) crystal structure 
with the room-temperature unit-cell parameters: a = 
15.2027 A, b = 15.1800 A and c = 13.5032 A (V = 
779.057 A 3 ) 3 The MA-K+ groups situated only in the 
large octahedral sites were found to be dynamically dis- 
ordered between eight equivalent orientations at high 
temperatures 32 Below the structural phase transition at 
T s = 220 K MA-K + groups become static and prob- 
ably ordered in an (anti)ferro-orientational order simi- 
lar to NH 3 K 3 C 60 .— Since the details of the low tem- 
perature structure are not known yet we base all our 
band-structure calculations on the experimental room- 
temperature structure. We select a single MA-K+ orien- 
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FIG. 1. (a) Band-structure and DOS for the room tempera- 
ture MAK 3 C 60 crystal structure. Red solid line is the interpo- 
lated band-structure using the maximally localized Wannier 
orbitals - see text for details. Black solid line: the total DOS; 
the red dot-dashed line, green dash-dot-dotted line, and blue 
dashed line represent the projected DOS on the first, second 
and the third Wannier orbital, respectively, (b) The split- 
ting of ti u bands at the V point, defining the orthorhombic 
crystal-field anisotropy. 



tation thus effectively imposing a ferro-oricntational or- 
der of the MA-K + groups. 

The three metallic t% u bands are well isolated from 
other molecular orbital-derived bands as is the case 
with every fullerene compound^ The band gap between 
"C 60 HOMO" (h u -derived bands) and "C 60 LUMO" (f lu - 
derived bands) is 0.96 eV, while the separation between 
"C 60 LUMO" and "C 60 LUMO+1" (ti g -derived bands) 
is 0.59 eV. For comparison we refer here to the corre- 
sponding gaps of 1.16 eV and 0.38 eV computed for cubic 
K 3 C 60 r^ and to the experimental values, roughly 1.8 eV 
and 1 eV respectively^ The underestimation of the band 
gaps, standard for LDA calculations, is fortunately of 
little consequence in our case. The neat separation of 
the different band groups allows us to focus entirely on 
ti u -derived bands close to the Fermi energy (Fig. [T]). 
Not surprisingly, bare DFT-LDA nonmagnetic calcula- 
tions yield for MAK 3 C 60 a metallic ground state with a 
half-filled £ lu band. However, several experiments proved 
that MAK 3 C 60 is an insulator ) 22 i 23 i 30 i 31 clearly suggest- 
ing that the mean-field DFT metallic state is driven to 
Mott-Hubbard insulating state by electron correlations. 
Since the phase diagram for non-cubic fullerides has al- 
ready been calculated within DMFT with included elec- 
tron correlations, it is sufficient at this stage to continue 
the characterization of MAK 3 C 60 in the non-correlated 
limit and then treat electron correlations, e.g. within the 
established MSDT phase diagram^ 

Orthorhombic crystal symmetry is reflected in the 
band-structure [Fig.QJa)] by (i) the removal of the three- 
fold degeneracy at the T point [Fig. HJb)]; and (ii) the 
inequivalence of Y and X points. At the T point we 



compute an energy splitting between the first and the 
second and between the second and the third ti u en- 
ergy levels of 5 = 71 meV and A = 140 meV, respec- 
tively. For comparison, the corresponding energy split- 
tings for NH 3 K 3 C 60 are S ps 30 meV and A ps 150 eV2& 
In general, A reflects the major orthorhombic crystal- 
field anisotropy and 5 tracks the smaller anisotropy in 
the ab plane. Surprisingly, band-structure anisotropics of 
the two systems are very much comparable despite some 
obvious crystallographic differences, like for instance the 
ratio c/a, which is 0.89 and 0.91 for MAK 3 C 60 and 
NH 3 K 3 C 60 , respectively. In order to position MAK 3 C 60 
on the MSDT fulleride (U, A) phase diagram we need 
to evaluate W, which is deduced directly from the com- 
puted band-structure. The resulting total t lu density of 
states (DOS), roughly shaped in three peaks, and shown 
in Fig.QJa) resembles that of NH 3 K 3 C 60 .— We note that 
the apparent disagreement between DOS's width and the 
band-structure splitting at the X point is due to very 
a small number of states at extreme values and it can 
thus be seen only in the zoomed view for small DOS 
values. We obtain a quite large DOS at the Fermi en- 
ergy, N(E-p) = 17 states/eV/Cg , which is a result of ex- 
panded lattice structure and the resulting smaller band- 
width W = 0.5 eV. For comparison, we refer here to the 
W ps 0.6 eV reported for NH 3 K 3 C 60 and K 3 C 60 , which 
are characterized by smaller unit cells^S 

The above results bring forward two interesting as- 
pects of the MA co-intercalation. The first is the po- 
sition of MAK 3 C 60 on the MSDT (U,A) phase diagram. 
Manini et. al£2- used A/W as a measure of the anisotropy, 
which also defines the distance to the MIT boundary in 
non-cubic fullerides. Based on A/W ps 0.25 NH 3 K 3 C 60 
is quite close to both 2- and 3- band metallic phases. It is 
thus expected that under pressure NH 3 K 3 C 60 would al- 
most instantly end up in one of these two metallic states. 
We find MAK 3 C 60 to be slightly more anisotropic, judg- 
ing from A/W = 0.29, thus deeper in the Mott-insulating 
phase and further away from the 3-band metallic phase 
(Fig.©. 

Application of a hydrostatic pressure, a standard ex- 
perimental method for increasing W and thus pushing 
systems across the metal-insulator boundary, has yet to 
be tried in this system to verify our calculations. How- 
ever, based on the above observation, we predict that 
larger pressures will be needed in MAK 3 C 60 to access 
metallic and possibly superconducting states. To esti- 
mate the required metal-insulator transition pressure in 
MAK 3 C 60 we compute the bandwidth as a function of 
decreasing unit-cell volume, mimicking the effect of an 
external pressure. The pressure dependence of the unit- 
cell parameters was taken from the high-resolution X-ray 
data measured under hydrostatic conditions^ To pre- 
vent unphysical contact between C 60 and MA groups we 
carried out a structural optimization for each volume. 
The bandwidth increases monotonically with decreasing 
unit-cell volume, or increasing pressure, with the slope 
of dW/dP = 100 mcV/GPa (inset Fig. |2J). If the effect 
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FIG. 2. MSDT phase diagram for fullerides adapted from 
Ref. I20I including MAK 3 C 60 position. Inset: Pressure depen- 
dence of the bandwidth computed for the room-temperature 
MAK 3 C 60 structure. The horizontal line marks the critical 
bandwidth 0.74 eV where (U/W) c = 1-35 and where metal- 
insulator transition is expected. 

of pressure could be reduced only to a pressure depen- 
dence of W, we would predict a steep increase of Tn in 
the insulating phase as a function of pressure, estimating 
for instance T N (1 GPa)/T N (0) - (W/W ) 2 = 1.44. At 
higher pressures insulator-metal instability is expected 
to occur. If we take critical ratio (U/W) c = 1.35 as 
appropriate for the orthorhombic structure^ and a typ- 
ical value for fullerenes U = 1 eV^ we predict that this 
transition should take place at around pmit = 2.3 GPa. 
Nonlinear effects on compressibility may push this criti- 
cal pressure slightly higher. 

Contrasted with the metal-insulator transition, mag- 
netic ordering cannot be understood purely in terms of 
overall properties such as the bandwidth or crystal-field 
anisotropy. In particular, the abnormally small Tn = 
11 K of MAK 3 C 60 as compared with T N = 40 K for 
NH 3 K 3 C 60 can hardly be justified by the fairly small dif- 
ferences either in the bandwidth or in the orthorhombic 
anisotropy. Accordingly, the lattice expansion alone can- 
not explain the small Tn : the antiferromagnetic ordering 
must be controlled by other degrees of freedom such as 
the interaction between co-intercaland molecules and C| 
anions and the related stabilization of a particular JT de- 
formation. We will test this hypothesis in the following 
sections. 



B. The Jahn- Teller effect 

Let us focus first on JT effect of Cg anion, which 
has been predicted 3 ^ but has been experimentally much 
more elusive in fullerene systems.— Orthorhombic crystal 
structures should provide a fertile ground for JT effect in- 
vestigations. In order to isolate JT effect we consider an 
artificial enlarged face centered orthorhombic Cg struc- 



ture using the room-temperature MAK 3 C 60 lattice pa- 
rameters multiplied by a factor of 1.5. To ensure charge 
neutrality, we add a uniform positive background in the 
DFT calculation. This positive background is not con- 
tributing to the orthorhombic crystal field, therefore any 
removal of t\ n degeneracy at the T point should arise 
solely from the JT effect on top of the weak residual crys- 
tal field from the periodically replicated Cg ions. The 
positions of C 60 carbon atoms were relaxed in order to 
obtain molecular distortions. The C|q ion deforms spon- 
taneously into a structure with Z?2h symmetry [Fig.|3Ja)]. 
This is the expected symmetry for t iu ® H g JT cou- 
pling involving H g vibrational modest - — Distortions 
are small with the maximum value of 2 pm. We estimate 
the energy scale for deformation by realizing that the re- 
laxed structure has AE t = 170 meV lower total energy 
than the starting structure with undistortcd icosahcdral 
C 60 . Former threefold ti u degeneracy of the LUMO is 
now removed, with ii u levels split equally by ~50 meV. 
The lowest £i u orbital is doubly occupied, while the third 
electron goes into the central ti u orbital, pinned at the 
Fermi level. The highest t± u orbital is empty. The to- 
tal energy difference, AE t , has several contributions: the 
JT effect, the crystal-field effect and bond- length correc- 
tion due to LDA approximation. To estimate only the 
JT energy scale we compare the above AE t with the one 
obtained for a structural relaxation with equal and fixed 
occupations of the ti u bands, which effectively hinders 
the JT effect. The difference in the total energy between 
these two calculations is 57 meV, a typical value for the 
JT cffect i 35 ' 37 Proper energy scale, typical size of defor- 
mations, the right symmetry of Cg the splitting and oc- 
cupation of ti u orbitals are strong indications that the 
observed distortion is indeed a result of the JT effect. 

The residual crystal field of C| ions is immediately 
seen for smaller lattice parameters, i.e. when room- 
temperature MAK 3 CgQ lattice parameters are for in- 
stance multiplied by a smaller factor of 1.25. Using the 
same procedure as above, the structural relaxation ended 
up with the same deformation of the Cgg molecule. The 
only difference is that the axis of JT deformation acci- 
dentally rotated from the crystal z- to y-direction. This 
indicates that we are dealing in both cases with the JT 
effect and that the change in the JT deformation axis is 
due to the existence of several equivalent minima of the 
lowest JT adiabatic potential energy surface (APES)* 3 " 5 - 
These minima become nearly degenerate for large lattice 
expansions. On the other hand, when lattice parameters 
are reduced down toward experimental MAK 3 Cg values 
the orthorhombic crystal field starts to play a role by 
making some of the APES minima deeper, thus promot- 
ing one specific JT deformation. 

Adding potassium atoms to the Cg structure results in 
an artificial orthorhombic K 3 Cg where even stronger or- 
thorhombic crystal field due to the close contact between 
the K+ and Cg ions are expected. The structural opti- 
mization of the C 60 carbon positions when starting from 
the JT distorted Cg atomic positions [Fig. [21a)] leads 
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FIG. 3. Radial molecular distortions of C50 ions in (a) or- 
thorhombic Cgp expanded structure, (b) K 3 C 60 orthorhombic 
environment. The size of blue circles and red squares is pro- 
portional to the amount of positive and negative deviations 
from the mean C 60 radius at a given carbon site, respectively. 

to deformations of the C 60 molecule shown in Fig. [HJb). 
Molecular deformations are slightly different because the 
crystal field additionally lifts the ti u degeneracy, hence 
producing pseudo-JT effect. Nevertheless, the resulting 
symmetry remains D211 and the maximal distortions of 
2.3 pm are similar to the previous cases. The same holds 
for the total-energy lowering, AE t = 140 meV. The JT 
effect is obviously still dominant over orthorhombic crys- 
tal field, which represents a smaller contribution to the 
total energy. The important message of this part is thus 
that the energy scale of the JT effect is Ejt ~ 60 meV 
and that for the experimental unit-cell volumes a compa- 
rably weak crystal field favors a particular Cg pseudo-JT 
deformation. 



C. Methylamine - C 60 interaction 

In the MAK 3 C 60 structure, methyl protons of the MA- 
K + groups approach Cg anions to very short distances 
and, based on this observation, suggestions about the 
weak hydrogen bond were formulated in the literature^ 
Such close contacts are expected to lead to rather strong 
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FIG. 4. C 60 radial deformations of MAK 3 C 60 for (a) room- 
temperature structure with relaxed C 60 atom positions only 
and (b) completely relaxed room-temperature structure. Blue 
spheres correspond to inward distortions of the C 60 cage and 
red boxes represent outward cage distortions. The sizes of the 
markers scale linearly with the amount of distortion. 



crystal fields and, according to the above discussion, to 
also affect strongly the pseudo-JT effect. It is not a priori 
clear which of the two effects is dominant in this struc- 
ture. Therefore, C|o deformations are now investigated 
in the room-temperature MAK 3 C 60 structure where only 
C 60 carbon atoms are allowed to relax. As expected, 
the presence of MA molecules has a dramatic effect on 
the Cg shape [Fig. H[a)]. The maximal cage distortions 
are significantly larger than in the previous cases - they 
reach up to 3.4 pm and the total energy is reduced by 
AE t = 303 meV during the structural optimization. In 
addition, even the I?2h symmetry of the distorted C 60 
molecule is lost. The maximal distortions are found for 
carbons facing methyl protons, at a closest-approach dis- 
tance of 226 pm. Dramatically larger distortions com- 
pared to those obtained for Cg or even for orthorhombic 
K 3 C 6 q structure indicate that the additional crystal field 
produced by MA-K + groups plays a dominant role over 
the JT effect in MAK 3 C 60 . 

We also address the possible hydrogen-bond formation 
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by relaxing all atomic positions, including those of MA- 
K+ groups. Hydrogen-bond traces can be detected on 
the tiny deformation of MA where C-H bond lengths of 
the CH 3 group are 110.3 pm, 110.4 pm and 110.6 pm, 
the last one corresponding to the hydrogen with the clos- 
est contact to the fullerene molecule. However, with the 
structural relaxation, the MA-K+ group rotates slightly 
away from the fullerene molecule, increasing the nearest 
H-C 6 q distance to 228.5 pm, i.e. by 2.7 pm longer than 
in the experimental structure [Fig. BJb)]. This distance 
still remains in the typical hydrogen-bond length range, 
but the Cgp - methyl proton contact elongation indicates 
that such bond must be very weak, if it exists at all. 

D. Magnetic interactions 

The above analyses point toward pseudo-JT effect and 
strong crystal fields in the MAK 3 C 60 arising from the 
MA-Cg interactions. The remaining question to be ad- 
dressed is how these effects influence the low-temperature 
magnetic properties. The appropriate starting point is 
to construct an effective tight-binding model for the t iu 
bands and then relate the transfer integrals to the ex- 
change coupling constants. To do so, we proceed by 
switching to the maximally-localized Wannier orbitals 
(WO)5 Setting the energy window to the range around 
the DFT t\ u band (Fig. [1]) we obtain three WOs, each 
one localized on the same C 60 molecule. Close to the 
fullerene cage, a typical WO (Fig. [5|) has the character- 
istic look of an appropriate combination of carbon 2p z 
orbitals^ The orbital of Fig. [5] resembles the DFT elec- 
tron density obtained at the T point. The distribution of 
2p 2 -like orbitals complies with the expected ii u symme- 
try. The other WOs (not shown) have a similar shape, 
but different orientations. Based on these WOs, we com- 
pute their on-site energies and transfer integrals between 
neighboring C 60 from the Kohn-Sham Hamiltonian as^ 

t aP {Ri) = (w(0,a)\H KS \w(Ri,l3)). (1) 

Here w(Ri, a) is a a-th WO and Ri is a lattice translation 
vector. For Ri = Eq. (TTJ) gives the on-site WO energies 
and for Ri ^ Eq. (fTJ) yields nine transfer integrals for 
each pair of neighboring C 60 molecules. By diagonalizing 
the tight-binding model, we verify that the obtained on- 
site energies and transfer integrals reproduce correctly 
the DFT electronic band structure as shown by red solid 
line in Fig. [Ha). 

We note that the computed on-site tight-binding 
Hamiltonian has non-zero off-diagonal elements, indicat- 
ing hybridization within the obtained set of WOs. By 
changing basis to an equivalent one where the on-site 
Hamiltonian becomes diagonal, we derive a new set of 
WOs with well-defined on-site energies. There are two in- 
dications that the new set of WOs is appropriate for our 
purposes: (i) the new WOs still reproduce band struc- 
ture and (ii) in the Mott-insulating limit such WOs are 




FIG. 5. A maximally localized Wannier orbital obtained for 
the MAK 3 C 60 structure. A isowavefunction surface at levels 
±65% of maximal value are depicted in red and blue color for 
positive and negative phase, respectively.— 



a good approximation for the localized ii u electronic or- 
bitals with well-defined on-site energies, in our case sep- 
arated by 32 and 67 meV. 

To estimate the exchange coupling constants we fo- 
cus only on the transfer integrals of the second WO, the 
half-filled one at intermediate energy. The neglect of the 
other (inter-band) hopping elements between neighbor- 
ing C 60 sites amounts to neglecting inter-band transi- 
tions. This selection can be justified by noting that (i) 
the retained WO is mainly responsible for the DOS at 
the Fermi level [see Fig. QJa)] and (ii) in a picture where 
the Coulomb repulsion is the dominating interaction, the 
intermediate orbital hosts the unpaired spin in the spin- 
1/2 state of each C|q molecular unit, with lower fully 
occupied and upper empty orbital. Within this approxi- 
mation we estimate the interfullercnc exchange coupling 
constants with Hubbard's expression Jy = At^/U , where 
we take U = 1 eV. For the considered MAK 3 C 60 struc- 
ture the strongest exchange couplings J3 = J4 = 1.9 meV 
are found for nearest C 60 neighbors along the (Oil) and 
(1 -1) directions. We stress that these J's are of the 
right order of magnitude as estimated from the measured 
Curie- Weiss temperature = 86 Ki2I All other exchange 
constants are much weaker, not exceeding 0.12 • J3, see 
Table H 

The obtained exchange network indicates a quasi-two- 
dimensional magnetic structure of MAK 3 C 60 . Low di- 
mensionality could account for the reduced experimen- 
tal Neel temperature. However, although this is cer- 
tainly a viable possibility, we stress that this conclusion 
is based on the high-temperature crystal structure. Since 
we proved that the interaction between MA and Cg is 
the governing factor for pseudo-JT effect, we stress at this 
point that the low-temperature exchange network may in 
fact be different. How sensitive interfullercnc exchange 
interactions are on the precise position of MA-K + groups 
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TABLE I. Exchange coupling constants for the experimental 
room-temperature structure [Fig. HJa)] and the relaxed struc- 
ture [Fig. 2£b)]. See text for details. 



Neighbors 


JExp. structure (meV) 


jRelaxed (meV) 


1: (1 1 0) 


0.03 


0.06 


2: (1 1) 


0.02 


0.02 


3: (0 1 1) 


1.94 


0.01 


4: (1 -1) 


1.89 


3.07 


5: (0 1 -1) 


0.00 


0.73 


6: (1 -1 0) 


0.22 


0.72 



becomes immediately evident when we consider relaxed 
DFT structure with a slightly rotated MA-K + groups 
[Fig. HJb)] . For this structure the strongest exchange is 
J4 = 3 meV, but now along (1 -1). The quasi low- 
dimensional picture seems to still hold since the next two 
strongest interactions are much weaker, i.e. J5 = 0.23 • J4 
(0 1 -1), and J 6 = 0.23 ■ J 4 (1 -1 0). All other exchange 
interactions are negligible (Table UJ). We stress that in 
all these structures the considered MA-K+ order is ferro- 
orientational. Any deviation from this MA-K + configura- 
tion, for example antiferro-orientational MA-K + order as 
the candidate for the low-temperature MAK 3 C 60 struc- 
ture, could significantly change the quasi low-dimensional 
character of our system. Moreover, if a certain degree 
of disorder in MA-K + orientations is present, then one 
would anticipate also a distribution of exchange coupling 
constants, and an even smaller ordered moment. That 
would explain the static magnetic order with a broad- 
local-field distribution measured by [iSR~2- This situ- 
ation is somehow reminiscent of ferromagnetic TDAE- 
C 60 , where disorder in JT deformations has been argued 
to be responsible for the reduced order parameter— ~— 
and broad distribution of local magnetic fields also mea- 
sured in /iSR experiment, for instance.— We conclude 
that pseudo-JT orientation is an essential factor in inter- 
fullerene exchange coupling in MAK 3 C 60 . Judging from 
the experimental high-temperature structure, MAK 3 C 60 
may even show low-dimensionality effects thus explain- 
ing the reduced Tn- However, other pseudo-JT orienta- 
tions may be active at low temperatures and even sup- 
pression of Tn due to the disorder effects cannot be ex- 



cluded at this stage. Further progress in understanding 
of magnetic properties of MAK 3 C 60 would be possible 
only when the low-temperature crystalline structure of 
MAK 3 C 60 is known. 



IV. CONCLUSIONS 

Antiferromagnetic Mott-Hubbard insulating 
MAK 3 C 60 had been studied by means of DFT cal- 
culations. We found the pronounced orthorhombic 
crystal field anisotropy in the electronic structure and 
judging from A/W = 0.29 MAK 3 C 60 should be deep in 
the insulating part of the phase diagram. Comparisons 
between different orthorhombic structures demonstrate 
that the JT effect is present and that its energy is around 
60 meV. However, the interaction between cointercalatcd 
MA molecules and Cg anions is comparatively strong 
in MAK 3 C 60 , leading to a strong crystal field and 
pseudo-JT effect stabilized at ~300 meV. The presence 
of strong pseudo-JT effect should be responsible for 
the suppression of Tn- Calculations based on the high- 
temperature structure with imposed ferro-orientational 
order of MA-K + groups suggest low-dimensionality in 
magnetic exchange; however other explanations such 
as a disorder in pseudo-JT cannot be excluded at this 
stage. A similar interplay of Mott and JT physics is 
believed to play an important role in the superconduc- 
tivity of underdoped cuprates. In the present study we 
emphasize the similar importance of JT effect in strongly 
correlated Cg structures and setting it as one of the 
most important parameters for the high-temperature 
superconductivity in fulleride salts. 
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